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ABSTRACT 

We show that the very close-by (19 pc) K0 star HD 189733, already found to be orbited by a 
transiting giant planet, is the primary of a double-star system, with the secondary being a mid- 
M dwarf with projected separation of about 216 AU from the primary. This conclusion is based 
on astrometry, proper motion and radial velocity measurements, spectral type determination and 
photometry. We also detect differential proper motion of the secondary. The data appear consistent 
with the secondary orbiting the primary in a clockwise orbit, lying nearly in the plane of the sky (that 
is, nearly perpendicular to the orbital plane of the transiting planet), and with period about 3200 
years. 

Subject headings: stars: low-mass, brown dwarfs - stars: individual: HD 189733 - stars: individual: 
HD 189733B - planetary systems - binaries (including multiple) 



1. INTRODUCTION 

Of the ~ 170 exoplanets in 146 planetary systems 
known at the present time 4 , the majority (~ 160) 
are revealed on ly by the reflex veloc i ty of their par- 
ent s tars, (e.g. iMavor fc Quelodll995t iMarcv fc Butlerl 
|1996|) . yielding their periods, eccentricities, semi-major 
axes, and msini values. However, a small number (9) 
of the known exoplanets transit their host stars, so that 
the inclination ambiguity is removed and - assuming a 
mass and radius for the primary - their actual mass and 
radius, hence their mean density may be determined. 

A number of the known exoplanets have been found 
to occur in multiple stellar systems. The small sam- 
ple of 22 such systems form an important class of ob- 
jects (for references, see: lEggenbereer et all2004bt f2005: 
Mugrauer et al.ll2004alhl.l2005alL Until now, however, no 
planet found in a multiple stellar system also transits its 
parent star. In this Letter, we note that the parent star 
of the recent ly-discovered transi ting extr a-sola r planet 
HD 189733b l)Bouchv et alJl2005l hereafter |B05^ is itself 
a member of a double-star system. 

In [JUwe present the evidence that the star HD 189733 
has a physical companion star. Specifically, "i'Z.W shows 
that it has a common proper-motion companion; i !2.2l 
shows that the companion's radial velocity is the same 
as HD 189733 within uncertainties; and i J2. 31 shows that 
the companion is a red dwarf star which must lie at ap- 
proximately the same distance as HD 189733. For these 
combined reasons the co-moving companion, now labeled 
HD 189733B, almost certainly must be a true physical 
companion. In Sj_D we detect a differential proper mo- 
tion, and obtain an initial estimate of the orbital motion 
of HD 189733B about HD 189733. Finally, discusses 
some implications of this finding, and avenues for future 
work. 
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2. EVIDENCE THAT HD 189733 HAS A PHYSICAL 
COMPANION STAR 

HD 189733 is a close-by (D = 19.3 ± 0.3pc) KO 
dwarf, with mas s of .82 ± 0.03M Q and othe r properties 
as de scribed in IB051 The 2MASS survey l)Cutri et alJ 
12001 lists a nearby red star, 2MASS 20004297+2242342 
(J = 10.12 ± 0.04, H = 9.55 ± 0.08, K s = 9.32 ± 0.03, 
J — Kg = 0.8), some 3.7 magnitudes fainter in Ks and 5 
magnitudes fainter in V (estimated from J, H, Kg). This 
star has an angular separation of 11.2" from HD 189733, 
lying 10.2" W and 4.6" to the S. 

HD 189733 is no t listed in the Wash ington Double Star 
Catalogue (WDS; IMason et afll200l() . but here we will 
show that 2MASS 20004297+2242342 is in fact its phys- 
ical companion; in anticipation of that result we hence- 
forth denote the companion as HD 189733B. Note that 
the latte r nam e is distinct from HD 189733b, the name 
given bv lB05l to the planetary companion to HD 189733. 
Our conclusion that the two stars form a bound binary 
system is based on the following evidence. 

2.1. Common Proper Motion 

HD 189733 has a relatively high proper motion of 
/i Q = -2.49 ± 0.68 masyr- 1 and us = -2 50.81 ± 0.53 
masyr _1 (Hipparcos; iPerrvman et al]|l997[) . To deter- 
mine whether HD 189733 and HD 189733B share com- 
mon proper motion, and hence may be members of the 
same physical system, we have inspected the follow- 
ing archival material: i) the digitized Palomar Observa- 
tory Sky Survey (POSS) scans 5 (POSS-I, 1951 R-band; 
POSS-II, 1990 R-band, 1992 B-band, and 1996 I-band); 
ii) the HST QuickV survey (1982, R -band), iii) the 
2MA SS 6 2000 J-, H- and K-band scans (|Skrutskie et alJ 
2000). In addition, in November 2005, we used the 
TopHAT telescope of the HAT Network at the Fred 
Lawrence Whipple Observatory (FLWO), to obtain I- 
band images. We also acquired 8 short exposure I-band 
frames in December 2005 with KeplerCam on the 1.2m 
telescope at FLWO. 

5 http:/ /archive. stsci.edu/cgi-bin/dss_form 

6 http://irsa.ipac.caltech.edu/ 
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Fig. 1. — HD 189733 in four different epochs on registered frames. The known proper motion to the South is well visible. The upper 
dashed line shows the declination in 1951, and the lower dashed line for 1999. The companion HD 189733B is not visible on the POSS-I 
frame (left panel) because HD 189733 is over-saturated and the scan resolution is not adequate. HD 189733B (indicated by B) is visible 
on the rest of the panels, displaced to the SW by ~ 11" from HD 189733; it is seen increasingly better from left to right, as the companion 
is relatively brighter at increasingly longer wavelengths (R, I and K bands). "X" denotes a faint star visible on POSS-I and the QuickV 
scans, before the Southward moving HD 189733B merges with it. "D" and "E" mark two faint stars that are visible on all of the frames, 
but one of the artificial filter glints ("G") merges with their position on the 2MASS frame. 



We used our home-grown fihat software environment 
(Pal et al. 2006, in preparation) to find sources on all 
the above images (~50 isolated, non-saturated stars), 
cross-match them by rejecting outliers, determine the as- 
trometric mapping between the images, and transform 
them to the same reference system. Visual inspection of 
the registered frames show i) the prominent Southward 
proper motion of HD 189733 (in harmony with the Hip- 
parcos values), ii) a much fainter co-moving companion 
that we identify as HD 189733B (for details, see Fig.QJ. 
The companion is clearly separated from HD 189733 on 
the 2MASS J, H and K scans (Fig. □ right panel), 
and is listed in the 2MASS point source catalogue as 
2MASS 20004297+2242342. Although not illustrated on 
Fig. n the co-movement is also demonstrated by the 
TopHAT and FLWO 1.2m I-band frames. No other co- 
moving companion is detected on these frames. 

To quantify the proper motion of HD 189733B, we have 
carried out astrometry on the QuickV, POSS-II, 2MASS, 
TopHAT and FLWO 1.2m observations. We used the 
2MASS catalogue as astrometric reference, where the 
quoted position uncertainty of bright, isolated sources 
is 120mas. We used the previously established pixel cen- 
ters of ~ 50 un-saturated stars that we found by fit- 
ting Gaussian profiles. By running our fihat/fistar 
star-finder algorithm on artificially generated frames, we 
found typical errors of the centroid positions to be on 
the order of 0.05pix (corresponding to 0.08" on QuickV, 
0.05" on POSS-II and 2MASS, and 0.1" on TopHAT and 
0.07" on FLWO 1.2m). We then derived the second or- 
der astrometric mappings between the X,Y coordinates 
and the 2MASS astrometric reference (HD 189733(B) 
were omitted from the fit), and used this to transform 
the pixel coordinates of the fra mes to the ICRS (a, 5; 
iSeidelmann fc Kovalevskvl2002|) system used by 2MASS. 
The rms around the fit was ~ 0.2" for QuickV, ~ 0.18" 
for POSS-II, - 0.05" for 2MASS, - 0.3" for TopHAT, 
and 0.1" for the FLW01.2m, respectively. 

The derived proper motion of HD 189733B is n a ,B = 
-4.1±9masyr -1 , /j, s>B = -264±12masyr _1 (see Fig.EJ, 
which are within 2a of the Hipparcos proper motion of 
HD 189733 itself. Therefore, it is clear that, within un- 
certainties, HD 189733 and HD 189733B share a com- 
mon proper motion and so are likely to comprise a bound 
system. 
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Fig. 2.— Proper motion of HD 189733B in RA (left) and Dec 
(right) relative to the 2MASS 1999 position. The two panels are 
on the same scale. The central solid lines show the linear fit to 
the data; the two other solid lines show the same fits using slope 
and intersection parameters differing by ±1(7. For reference, the 
dashed line shows the relative proper motion of HD 189733 from 
Hipparcos, offset by 1.0 arcsec for clarity. 

2.2. Radial Velocity 

A common radial velocity is a further indication that 
the two stars are physically associated. In order to test 
this possibility, we obtained spectroscopic observations of 
both the primary star and the suspected companion, us- 
ing the Center for Astr ophysics (CfA) Digital Speedome- 
ter (DS; lLathamll992[) at the 1.5m Tillinghast telescope 
of FLWO, Arizona. 

Seven DS observations have been made of the star 
HD 189733 dating back to 1995, with the two most re- 
cent being 2005 Dec 10 and Dec 17. For each of these 
a radial velocity was o btained on the CfA Native Sys- 
tem velocity reference ijStefanik et al.lll999j) . The mean 
and standard deviation of these measurements is V^ a( j 
= —2.38 ± 0.20 kms -1 , with no significant evidence for 
a long-term velocity variation over the past 10 years. 
For HD 189733B, two DS observations were made, on 
2005 Dec 10 and 17, yielding a mean radial velocity 
^rad.B = —3.1 ± 1.0 kms -1 . Thus, within observational 
uncertainties the difference in measured velocities is con- 
sistent with the two stars being physical companions. 



A stellar companion in the HD 189733 system 
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2.3. Characteristics of the Star HD 189733B 

There is still a remote chance that the apparent com- 
panion could be a distant giant star with high tangential 
velocity or a very close-by low luminosity and velocity 
sub-dwarf, that happens to have the same sky position, 
radial velocity, and proper motion. 

While the DS has the primary goal of radial-velocity 
measurements, corr elation of the sp ectrum with template 
spectra based on the lKuruc3 l)1993(l models also yields in- 
formation on T c ff, \ogg and stellar rotation usini. We 
obtained logg = 4.5 ± 0.3 - a value appropriate to a 
main sequence M dwarf. Correlation with observed spec- 
tra of M dwarfs with spectral type ranging from MO.O to 
M5.5 gives best correlation near M3.5. We also observed 
HD 189733B with the FAST spectrograph on the FLWO 
1.5m telescope. Visual comparison of the resulting spec- 
tra to known M dwarf spectra suggests a spectral type of 
M4V. Based on its spectral type, plus apparent magni- 
tude, we then infer a distance to HD 189733B consistent 
with the 19.3 pc distance to HD 189733 . 
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Fig. 3.— Location of HD 189733B on the IBaraffe et alJITflflSl) 
5Gyr isochrone with other M-d warf binary companions o f stars 
with known planets, plotted from Eeecnbcrecr ct al. (2004b). Stel- 
lar masses are labeled along the isochrone in units of Mq. 

We independently estimated the spectral type of 
HD 189733B from 2MASS photometry. Because of 
the slightly overlapping profile of the close-by, bright 
HD 189733 we performed aperture photometry of 
HD 189733B on the 1999 2MASS scans after subtract- 
ing off the Gaussian profile of HD 189733, and using 
~ 50 isolated stars with original 2MASS photometry 
as reference. This analysis yields J = 10.147 ± 0.02, 
H = 9.551 ± 0.03, and K s = 9.318 ± 0.02. These val- 
ues are within 0.03mag of the 2MASS Point Source Cat- 
alogue values, but we find smaller errors, and slightly 
redder J — K color index. We transformed our mea- 
sured J — Ks value (namely 0.829 ± 0.03), t o J — K on 
the B essell and Brett (BB) system following iCarpenterl 
{2 QQ 5jl ; to obtain (J - K) BB = 0.864 ±0.04. Then we de- 
rived the absolute magnitude Mh (8.13 ± 0.045) from 
our measured H magnitude and an assumed distance 
equal to that of HD 189733 an d plotted these values o n 
the color-magnitude diagram of M ugrauer et all (12 005a V 
which also plots a 5-Gyr isochrone from Baraffc ct al. 
(1998). The position on the color magnitude diagram 



(Fig. |3J| corresponds to a s tellar mass of 0.175 to O.2M , 
which according to iCoxl ((2QQQ/) corresponds to an M 
dwarf with spectral type of about M5. The good fit to 
the isochrone supports our assumption that the distance 
to HD 189733B is similar to that to HD 189733. 

From the above analyses we conclude that 
HD 189733B is an M dwarf with spectral type in 
the range M3.5 to M5, with common proper motion, 
radial velocity, and distance to HD 189733. The relative 
positions (from astrometry, and assuming equidistance), 
and relative velocities of the two stars (from proper 
motion and radial velocity data), along with the IMq 
total system mass, are consistent within 2a with the 
system being gravitationally bound. Although it would 
formally be possible for an interloper M dwarf star to be 
passing very close to HD 189733 at the current epoch, 
with a relative speed so small as to make it almost 
gravitationally bound, this possibility is so remote that 
we conclude that the two stars do indeed form a bound 
system with projected separation about 216 AU. 

3. ORBITAL CHARACTERISTICS OF THE SECONDARY 

If the true separation of HD 189733 and HD 189733B 
is close to the projected separation, and the orbit is cir- 
cular, and the to tal sy stem mass is ~1M Q (O.82M for 
HD 189733 from IBOsT and O.2M for HD 189733B from 
H2.3(l then the orbital period is ~3200 years, correspond- 
ing to ~2kms -1 orbital motion. For a face on orbit this 
would yield an observable 22 masyr -1 differential proper 
motion in addition to the co-movement. 

In an attempt to detect this, we used two pairs of 
2MASS images in the H and K bands obtained in 1999 
and 2000 (one of these is the right-hand image in Fig.^l, 
with mean epoch 2000.073, and compared them with the 
8 high resolution I-band frames of the field taken with 
the FLWO 1.2m telescope at epoch 2005.970. We note 
that although the other data were useful in confirmation 
of the common proper motion (as shown in H2.1fl . because 
of the saturated HD 189733 image (POSS, Quick- V) and 
low S/N (TopHAT) they were not used in this precise as- 
trometry aimed at refining the differential proper motion. 
The astrometry of HD 189733 over the 5.897 year base- 
line yields a proper motion of [i a = —8 ± 5 masyr -1 and 
lis = —245 ± 8 masyr" 1 . Because this agrees with the 
more precise Hipparcos value within uncertainties ft |2.1f> . 
we adopt that value for the proper motion of HD 189733. 
For HD 189733B we find fi a ,B = —3 ± 5 masyr -1 and 
[isb = —272 ± 5 masyr -1 , which imply a differential 
proper motion relative to HD 189733 of A/i Q = — 1 ± 5 
masyr -1 and A/x,5 = —21.2 ± 5 mas yr -1 (Fig. 

Formally, the data indicate a detection of relative 
proper motion at the 4<t level. The position of the com- 
panion, and its direction and magnitude of relative mo- 
tion, arc consistent with orbital motion in a clockwise 
orbit roughly in the plane of the sky. 

However, there could be underlying systematic effects 
due to the different instruments and bandpasses used 
for the earlier-epoch 2MASS data and the later-epoch 
FLWOl.2 data. We estimated one of these, namely 
the effect of stellar profile merging of HD 189733B 
with HD 189733 that is different on the 2MASS and 
FLWO 1.2m frames. By subtracting off the Gaussian pro- 
file of the primary, we found that the derived proper mo- 
tion of HD 189733B changed by only 2masyr -1 . Nev- 
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ertheless, while the detection of orbital motion roughly 
in the plane of the sky seems secure, it is premature to 
derive specific orbital parameters for the system. 



-4.70 



-4.75 



-4.80 - 



2 -4.85 

Q 
X 

e 

o 

S -4.90 
> 



-4.95 
+10.00 







□ . 


' 2 MASS D 






D f 








□ 






/ 


FLW01 .2 






/ ° 








A 1 

/o '8? 










o °o 











+10.05 +10.10 +10.15 +10.20 
DX (from HD1 89733, West, in arcsec) 



-10.25 



Fig. 4. — Differential proper motion of HD 189733B relative to 
HD 189733. Open squares and circles indicate individual 2MASS 
(epoch 2000.073) and FLWOl.2 (epoch 2005.970) data points re- 
spectively. Filled square and circle represent the mean of those 
respective data points, with error bars representing 1 standard 
deviation of those means. Dashed line depicts a circular orbital 
path about HD 189733 in the plane of the sky, passing through 
the 2MASS position. The triangle depicts the position at epoch 
2005.973 for such an orbit. 



4. DISCUSSION 

Several studies have been done on the characteristics 
of close-in planets orbit ing the primary st a r of a mul- 
tiple star system. Thus lEggenberger et all l)2004ri and 
references therein) found a tendency for massive plan- 
ets (msini greater than 2Mj) to occur preferentially in 
multiple-star systems. However, HD 189733b is a low- 
mass planet (1.1 5 Mj without sini ambiguity), and yet it 
is in a multiple stellar system, thus weakening the dis- 
tinction between single and multiple stars. 

HD 189733b exceeds all other known extrasolar plan- 
ets in binary systems in its proximity to its parent 
star (a p ;=0.031 AU). However, a number of planets in 
multiple-star systems are nearly as close (e.g. Tau Boo 
b, a p; =0.05 AU; HD 75289b, a p; =0.046 AU). 

The data suggest (see S£$J) that the binary orbit is 



likely to be nearly face on, i.e. the orbital plane would 
be nearly orthogonal to the orbital plane of the planet 
HD 189733b, which by virtue of its transit we know to 
have an inclination of nearly 90°. A detailed calcu- 
lation based on the relative velocities and positions of 
the two stars shows that the orbit of HD 189733B and 
HD 189733b cannot be coplanar at the 4-cr level (Pal 
et al, in preparation). When additional transiting plan- 
ets are discovered in multiple star systems, it should be 
possible to study the relation of the system architecture 
(e.g. mass, semi-major axis of the stellar secondary) to 
planet properties in better detail. 

Because the HD 189733 system is so close-by, it should 
be possible to do excellent astrometry over the next few 
years with modern high precision astrometric techniques 
from ground or space. The changing radial velocity signal 
of either or both of HD 189733 and HD 189733B might 
also be detectable after a few years of monitoring with 
current high-precision radial velocity devices. High pre- 
cision proper motion and radial velocity measurements 
should lead to a good characterization of the orbital char- 
acteristics of the binary system and their relation to the 
orbit of the transiting planet. 

Finally, we note that because this system is only 19 
pc from Earth and hence both stellar components are 
unusually bright for their spectral types, many addi- 
tional follow-on observations requiring high resolution 
spectroscopy or high precision photometry will be fea- 
sible. This should permit a full characterization of the 
system including the possible detection of additional low 
mass components. 



This work was funded by NASA grant NNG04GN74G. 
Work by G. A. B. was supported by NASA through grant 
HST-HF-01170.01-A Hubble Fellowship. A. P. wishes to 
acknowledge the hospitality of the Harvard-Smithsonian 
Center for Astrophysics, where part of this work has been 
carried out. Work of A. P. was also supported by Hungar- 
ian OTKA grant T-038437. D. W. L. and R. P. S thank 
the Kepler mission for support through NASA Coopera- 
tive Agreement NCC2-1390. We acknowledge the use of 
the 2MASS survey frames, and the Palomar Sky Survey 
digital scans. We thank Perry Berlind at FLWO for tak- 
ing the FAST spectra of HD 189733B, and Emilio Falco 
for the follow-up exposures using the 1.2m telescope. 



REFERENCES 



A. N. Cox. 2000, Allen's astrophysical quantities, 4th ed. Publisher: 

New York: AIP Press; Springer, 2000 
Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, 

A&A, 337, 403 
Bessell, M.S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 
Bouchy, F., et al. 2005, A&A, submitted, astro-ph/0510119 
Carpenter, J. 2005, ApJ, 121, 2851 
Charbonneau, D. et al. 2005, astro-ph/0508051 
Cutri, R. M., et al. 2003, VizieH Online Data Catalog, 2246, 
Eggenberger, A., Udry, S., Mayor, M., Beuzit, J.-L., Lagrange, 

A. M., & Chauvin, G. 2004, ASP Conf. Ser. 321: Extrasolar 

Planets: Today and Tomorrow, 321, 93 
Eggenberger, A., Udry, S., & Mayor, M. 2004, A&A, 417, 353 
Eggenberger, A., Mayor, M., Naef, D., Pepe, F., Queloz, D., Santos, 

N. C, Udry, S., Lovis, C. astro-ph/0510561 
Kervella, P., Thevenin, Di Folco, K., & Segransan, D. 2004, A&A, 

426, 297 

Kurucz, R. 1993, ATLAS9 CD-ROM No. 13. Cambridge, MA, 
SAO, 1993., 13 



Latham, D. W. 1992, ASP Conf. Ser. 32: IAU Colloq. 135, 32, 110 

Marcy, G. W., & Butler, R. P. 1996, ApJ, 464, L147 

Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & 

Worley, C. E. 2001, AJ, 122, 3466 
Mayor, M., & Queloz, D. 1995, Nature, 378, 355 
Mugrauer, M., Neuhauser, R., Seifahrt, A., Mazeh, T., & Guenther, 

E. 2005, A&A, 440, 1051 
Mugrauer, M., & Neuhauser, R. 2005, MNRAS, 361, L15 
Mugrauer, M., Neuhauser, R., Mazeh, T., Guenther, E., & 

Fernandez, M. 2004, Astronomische Nachrichten, 325, 718 
Mugrauer, M., Neuhauser, R., Mazeh, T., Alves, J., & Guenther, 

E. 2004, A&A, 425, 249 
Perryman, M. A. C, et al. 1997, A&A, 323, L29 
Seidelmann, P. K., & Kovalevsky, J. 2002, A&A, 392, 341 
Skrutskie, M. F., et al. 2000, VizieR Online Data Catalog, 1, 2003 
Stefanik, R. P., Latham, D. W., & Torres, G. 1999, ASP 

Conf. Ser. 185: IAU Colloq. 170, 185, 354 



